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Abstract 

Point multipole expansions are widely used to gain physical insight 
into complex distributions of charges and to reduce the cost of computing 
interactions between such distributions. However, practical applications 
that typically retain only a few leading terms may suffer from unacceptable 
loss of accuracy in the near-field. We propose an alternative approach for 
approximating electrostatic charge distributions, Optimal Physical Multi- 
poles (OPMs), which optimally represent the original charge distribution 
with a set of point charges. By construction, approximation of electro- 
static potential based on OPMs retains many of the useful properties of 
the corresponding point multipole expansion, including the same asymp- 
totic behavior of the approximate potential for a given multipole order. 
At the same time, OPMs can be significantly more accurate in the near 
field: up to 5 times more accurate for some of the charge distributions 
tested here which are relevant to biomolecular modeling. Unlike point 
multipoles, for point charge distributions the OPM always converges to 
the original point charge distribution at finite order. Furthermore, OPMs 
may be more computationally efficient and easier to implement into exist- 
ing molecular simulations software packages than approximation schemes 
based on point multipoles. In addition to providing a general framework 
for computing OPMs to any order, closed-form expressions for the lowest 
order OPMs (monopole and dipole) are derived. Thus, for some practi- 
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cal applications Optimal Physical Multipoles may represent a preferable 
alternative to point multipoles. 

1 Introduction 

Point multipole expansions are widely used to gain physical insight by providing 
a simplified expression for a complex distribution of sources of potential fields, 
such as electrostatic potential due to a charge distribution. Point multipoles 
provide a means of decoupling the underlying features of a source distribution 
from the observation point. Thus, one can obtain physically meaningful insights 
into the macroscopic properties of the distribution, such as the familiar dipole 
moment. 

In addition to having theoretical utility, the point multipole expansion has 
been used to simplify practical calculations. For example, algorithms such as the 
fast multipole method [TT], use point multipoles to reduce the computational 
complexity of calculating pairwise interactions between large charge distribu- 
tions. Let Rq be the distance of the furthest charge in the distribution from 
the expansion center. Since each successive term in the multipole expansion 
decays more rapidly with distance than the previous term, at large distances 
R ^ Rq the series converges quickly and so these methods are able to obtain 
reasonably accurate results by keeping only the lowest order terms. However, 
at distances not much larger than Rq, the accuracy of the approximation dete- 
riorates quickly. Since, in practice, the potential often needs to be calculated in 
regions where the assumption R^ Rq does not hold (the near- field), the point 
multipole expansion is suboptimal for many practical calculations. For exam- 
ple, in atomistic molecular simulations, amino acids interacting inside a single 
protein are rarely more than several times Rq apart. For these calculations, 
point multipoles can be expected to provide a suboptimal approximation of the 
original distribution. 

We investigate an alternative to the point multipole expansion for approx- 
imating charge distributions - Optimal Physical Multipoles. Optimal Physical 
Multipoles approximate a charge distribution by a small number of point charges 
spatially distributed in an "optimal" way, to be precisely defined below. Since 
OPMs have a finite size, they may provide better representation of the original 
spatially extended charge distribution, i.e. a more accurate representation of 
the potential in the near-field. Consider for example the extreme case where the 
charge distribution consists of 6 charges of equal magnitude spaced equidistant 
on a line and alternating negative to positive (Figure [1]). For such a distri- 
bution, an optimally placed point dipole still produces more than 25 times the 
average error produced by an Optimal Physical Dipole at a distance of 2Rq from 
the center of geometry of the original charge distribution. Furthermore, later 
we will prove that for at least the lowest order OPMs, i.e. Optimal Physical 
Monopole and Optimal Physical Dipole, OPMs are at least as accurate as an 
equivalent order point multipole. From a practical standpoint. Optimal Physical 
Multipoles may be easier to implement in applications that already utilize point 
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Figure 1: Relative errors in electrostatic potential around an example net 
neutral charge distribution consisting of 6 charges with equal magnitude 
spaced equidistant on a line and alternating negative to positive, see inset. 
An Optimal Physical Dipole produces less than 3.5% error for all angles 
(solid purple line), whereas an optimally placed point dipole produces over 
45% error at some angles (dashed green line). The Optimal Physical Dipole 
consists of a negative and positive charge, of equal magnitude, placed on 
the left and right ends of the original charge distribution respectively. On 
the other hand, an optimally placed point dipole, as defined in ^ 117] . 
is located at the center of geometry. The relative error is computed as 

(^|$(2i?o, 0) - $.e/(2i?o, 0)\ /^Jil/An) Jo " /;($„/(2i?o, e))2 sm{0) d(f,d0)^ x 

100% where $(2i?o, ^), *&re/(2i?o, ^^) are the electrostatic potentials of the 
approximate and the reference (original) charge distributions respectively. 
The electrostatic potentials are computed at a distance 2Ro from the center 
of geometry, where 2Rq is the size of the charge distribution (twice the 
distance from the center of geometry to the outermost charge in the original 
distribution). 



charges, i.e. in many molecular dynamics packages [161 118] and are potentially 
more computationally efficient than point multipoles. 

Representing charge distributions by a small number of point charges is not 
by itself a novel idea. There are a number of methods, such as RESP [4], 
CHELP [6], CHELPG 'F, CHELMO [l9], coarse graining O [El and others 
[23] that empirically fit a set of point charges to a given charge distribution 
by minimizing various error metrics in electrostatic potential over some volume 
or surface surrounding the charge distribution. A key difference between these 
methods and Optimal Physical Multipoles introduced here is that OPMs inherit 
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the physically appealing asymptotic properties of the point multipole approxi- 
mation, i.e. the error in potential falls off at least as fast as l/it!'^+^, where R is 
the distance from the origin and k is the highest order of the multipole terms 
retained in the expansion. 

The rest of this work is organized as follows. We first review the multipole 
expansion concept. Next, we describe the theoretical basis for Optimal Physical 
Multipolcs. We then use this theoretical formalism to derive closed-form expres- 
sions for the Optimal Physical Monopole and Dipole, discuss their properties 
and compare their accuracy with the corresponding point multipoles. Finally 
we provide a cost-benefit analysis of Optimal Physical Multipolcs versus the 
point multipole method and discuss the prospect of expansions to higher or- 
ders. Potential uses and future work is discussed in "Conclusions". Although 
point multipolcs have been used to study a wide range of potentials, this proof- 
of-concept study focuses on just one of the most common application of the 
concept - electrostatics. Nevertheless, we believe that the ideas presented here 
will be generally applicable to other 1/R, Coulomb- like, potentials as well. 



Here we will give a brief overview of the familiar formalism of the point multi- 
pole expansion. Since many practical applications, such as molecular dynamics 
simulations, use point charges, for notational simplicity we will consider discrete 
charge distributions, but the results also hold for continuous distributions. 

Consider a set of N point charges Qn {n = 1, 2, iV) located at positions r„ 
around some chosen origin. Then the potential ^{R), of this distribution at a 
point R from that origin is given by the familiar Coulomb Potential 



For distances R > Rq where R = \\R\\ and i?o = waxdlr^H), a Taylor se- 
ries expansion of the potential above gives the classic multipole expansion. In 
Cartesian coordinates we obtain 



2 Multipole Expansion 
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and q,p,Q,0 are known as the monopole, dipole, quadrupole and octupole 
moments respectively, and Sij is the Kronecker delta. The multipole moments 
are symmetric tensors where the lowest order non-vanishing multipole is origin 
independent. 



3 Definition of Optimal Physical Multipoles 



For a given set of original charges qi {i = 1,2, N), we want to determine 
a smaller representative set of charges qk {k = 1,2, ...,K < N) such that the 
potential due to these representative charges, ^{R) best approximates the po- 
tential of the original distribution, $(-R). In the following sections, we will 
outline the general method for determining this optimal set of representative 
charges, i.e. the Optimal Physical Multipole representation. 

3.1 The Error Metric 

Determining the best representative charge distribution is critically contingent 
upon the definition of the error metric used. In general we are concerned with 
obtaining the best representation of the original potential at any arbitrary point 
in space outside the distribution. Thus, for the error metric. A, one typically 
chooses the root mean square (RMS) of the error in potential over some volume 
V excluding the volume Vq containing the charge distribution being approxi- 
mated, i.e. 



In principle, one can derive the optimal charge placement fn} by minimiz- 
ing the integral given in Eq. ^ with respect to the values of the new charges, 
{q-n} and their positions However, as the number of charges in the rep- 

resentative distribution grows, this equation can be difficult to minimize, let 
alone to find closed-form, analytic expressions for the placement and magni- 
tude of the charges composing the representative distribution. Furthermore, 
charges chosen in this manner are not guaranteed to have the same multipole 
moments as the original distribution [19] . This can lead to misinterpretation of 
the properties of the distribution and, potentially, to unphysical results. At the 
very least, we would like the new representation to inherit the same transparent 
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asymptotic behavior of the corresponding point multipole expansion, but with 
greater accuracy expected from an extended distribution that can better mimic 
the original. To simphfy the problem, we recast Eq. ([7]) in spherical coordinates 
and consider the error inside a spherical shell centered on the chosen multipole 
expansion center, and with arbitrary outer radius R > Rq, where Rq is defined 
as before, i.e. the distance from the expansion center to the outermost point 
charge. The error metric now becomes 



= ^^(j^3 _ / / / - HR)\ R'^HmmdR (8) 



i?o 



where and (j) are the usual spherical coordinate inclination and azimuth angles. 
In spherical coordinates, the multipole expansion is given by 

^' e=om=-e 

where Y™ are the standard spherical harmonics, * is the complex conjugate, 
are the spherical multipole moments, and £ is the multipole order. 

N 

qT^J2l^^nyr{0n,cl^nr (10) 
n=l 

Using this expansion as our error metric, Eq. ([8]), becomes 

-Z. 27TT ' 
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where and are the spherical moments of the original and representa- 
tive charge distributions respectively. Since the spherical harmonic integral 
!T \yrid,(l))\^sm{8)ded(j} = 1, Eq. (dll), can be further simplified to the 
following form |17) , 

q ^ OO m—i 

A' = i / y = y \qT - Wl'dR (12) 

Ha 

From the structure of Eq. (jl2p we see that minimizing the difference between 
the successive multipole moments of the original and Optimal Physical Multipole 
distributions is equivalent to minimizing the total error, i.e. by minimizing 
each term in the error expansion we minimize the RMS error in electrostatic 
potential. Note that the procedure does not depend on the parameter R, and 
thus the method does not require explicit integration over a given region. The 
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use of the multipole expansion in this way allows for the sought after distinct 
separation of terms by the rate at which they decrease as a function of R, i.e. 
the monopole term falls off as the dipole falls off as etc. A representation 
which makes terms up to order k in Eq. (I12p equal to zero will produce total 
error whose leading term falls off as 

3.2 Optimal Physical Multipole Approximation 

We define an Optimal Physical Multipole of order /c to be a set of charges 
satisfying the following conditions: 

1 . It is the smallest set of point charges which make every term in the error 
expansion, Eq. (|12p. through order k equal to zero (for the given charge 
distribution being approximated). 

2. If the charge magnitudes and positions are not uniquely determined by 
condition 1, the charges are chosen such that the k + l,k + 2,... terms 
in the error expansion are sequentially globally minimized until all the 
charges are uniquely defined. 

Optimal Physical Multipoles as defined above have the useful property that 
for any point charge distribution they will, at some order fc, reproduce the 
original point charge distribution. For, at some order fc, the 0PM will require 
the same number of charges N as the original distribution to satisfy the first 
condition of the Optimal Physical Multipole definition. Furthermore, the 0PM 
must be made identical to the original distribution to satisfy condition 2 from the 
definition above. To see how this arises note that the error expansion contains 
terms of the form (q™ — g™)^- Thus, an Optimal Physical Multipole which is 
identical to the original distribution will have q™ = for all £, m and clearly 
having = for all m will be a global minimum of the error expansion 
term of order £. Since making the 0PM identical to the original distribution 
will always provide a global minimum to all orders, setting the 0PM equal 
to the original distribution provides the only way to place the charges such 
that they satisfy the second condition defining an 0PM. We note that the 
convergence of Optimal Physical Multipoles to the original distribution after a 
finite number of computations sets OPMs apart from point multipole expansions 
which generally require an infinite number of terms to exactly reproduce a given 
charge distribution. 

The minimizations of the error metric in Eq. ()12[) which are required to 
define an Optimal Physical Multipole can be done numerically to arbitrary or- 
der. This numeric procedure for calculating the Optimal Physical Multipole 
representation may be particularly useful in situations where the charge distri- 
butions are relatively static and thus the optimal representation does not need 
to be recalculated at each time step of a given simulation. Ideally, however, one 
would like analytic expressions that can be used to compute OPMs at a reduced 
computational cost. 
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4 Analytic Expressions for Low Order Optimal 
Physical Multipoles 



In the following sections, we apply the general methodology developed in section 
3, to derive low order OPMs and test their accuracy. In sections 4.1 and 4.2, 
we obtain analytic forms for the lowest order Optimal Physical Multipoles in 
two important special cases, namely, we define the Optimal Physical Monopole 
for systems with non-zero charge and the Optimal Physical Dipole for systems 
with zero charge but non-zero dipole moment. These analytic expressions not 
only provide physical insight but are more computationally efhcient than the 
numerical minimization procedures that are in general required to obtain an 
approximate charge distribution. Thus, these analytic expressions may be par- 
ticularly useful in applications such as molecular dynamics where computational 
speed is critical. 

4.1 Optimal Physical Monopole 

The Optimal Physical Monopole, order fc = 0, consists of a single charge. As 
long as the single charge has magnitude q = q, i.e. is equal to the total charge of 
the original distribution, the monopole term of the error expansion will obviously 
be zero. Thus, any charge with magnitude equal to the net charge of the original 
distribution will satisfy the first condition which defines a Optimal Physical 
Multipole. Now, the remaining parameters, namely the position of the charge, 
are chosen to satisfy condition 2 of the Optimal Physical Multipole definition, 
namely to minimize the dipole term in the error expansion. In this particular 
case, the A: -I- 1 term of the error, i.e. the dipole term, can be made identically 
zero by solving 



for X, y, z where Px,Py,Pz are the x, y, z components of the dipole moment p of 
the original distribution. Solving the above equations we have 



So, a charge of magnitude q placed at f satisfies both conditions which define 
Optimal Physical Multipole of order fc = 0. 

The Optimal Physical Monopole can offer substantial improvements com- 
pared to the use of a single charge placed at the common choice of the geomet- 
ric center. In particular, there are biologically relevant examples, see Fig. 
where a point monopole at the origin produces an RMS error 5 times greater 
than that of the proposed Optimal Physical Monopole. 



Px - q- X = 

Py-q-y = 

Pz-q - z = 



(13) 
(14) 
(15) 



P 

r — — 

q 



q = q 
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Figure 2: Relative error in electrostatic potential produced by an Optimal 
Physical Monopole representing the charge distribution of an amino acid that 
carries a net charge of +l|e| (glutamic acid at physiological pH). The inset image 
shows the plane bisecting the molecule which defines the angle plotted in the 
figure; to illustrate the original charge distribution, the electrostatic potential 
of the glutamic acid is plotted on its molecular surface as a color map with blue 
representing positive and red representing negative potential |8]. The relative 
error is computed as in Fig. 1, (ratio of the absolute error to the RMS average of 
the reference potential), at a distance 2Rq which as before is twice the distance 
from the center of geometry to the outermost charge. The Optimal Physical 
Monopole, located at the center of charge, produces less than 11% error (solid 
purple line) for all angles whereas the point monopole, located at the center of 
geometry, produces over 45% error (dashed green line) for some angles. 



Despite its similarity to the center of mass, the existence of negative charge 
makes the center of charge term fundamentally different from center of mass. In 
particular, as the total charge in the distribution approaches zero, the center of 
charge tends towards infinity. However, in practice, charge is discrete, so there 
is a limit to the maximum distance the center of charge can fall outside of a 
given distribution. 

4.2 Optimal Physical Dipole 

For an uncharged distribution, the Optimal Physical Dipole consists of two 
charges Qi — q and q2 = —Q located at positions ri and r2 respectively. Thus, 
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it takes 7 parameters to uniquely define an Optimal Physical Dipole. By setting 

N 

q{fl - ^2) = ^ qnVn (18) 

n=l 

the dipole term in the error is zero and the two charges satisfy condition 1 of 
the Optimal Physical Dipole definition. Now, we will rewrite the positions fi 
and f2 in the following form: 

ri = d + — —- 

2q 

r, = d-^— (19) 

where d represents the geometric center between the two charges. We can see 
that these positions satisfy relation (|18p automatically. By writing the positions 
of the charges in this manner, we have divided the process of determining the 
remaining parameters which define the Optimal Physical Dipole into two steps, 
namely, finding the optimal placement of the charges, d, and finding the optimal 
magnitude of the charge, q. Note that finding the optimal charge value fixes the 
separation between the two charges, since the dipole moment of the represen- 
tative distribution has been constrained to equal the original dipole moment. 



The placement of the geometric center of the charges composing the Optimal 
Physical Dipole, d, that minimizes the quadrupole term of the error expansion 
is given by 




'^i,j=x,y,z QjiPiPj 
4p2 



Pk (20) 



where k = x,y,z and thus the d^'s are the components of d. This optimal 
position, known as the center of dipole, was derived previously [17] for a differ- 
ent purpose, namely for matching point multipole expansions between different 
charge distributions. Now, unlike the point dipole, the Optimal Physical Dipole 
has physical size and thus an additional parameter with which to further min- 
imize the error with respect to the given potential. In other words, Eqs. (|20p 
and (|18p . determine only 6 of the 7 parameters required to define an Optimal 
Physical Dipole. Since the quadrupole moment is the lowest order non-zero 
term remaining in the error expansion, by choosing the optimal charge value, 
we want to further minimize the quadrupole term in the error. However, for 
any charge value q an Optimal Physical Dipole placed at the center of dipole 
has no quadrupole moment as can be seen by setting N = 2, substituting the 
center of dipole, Eq. ([20]) and qi = —q2 = q into Eq. (fT9|) then substituting 
these variables into Eq. ([6]). Thus, the quadrupole term in the error, Eq. (fT2|) 
is unaffected by the choice of the charge magnitude, q and the quadrupole term 
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Figure 3: Relative error in electrostatic potential produced by Optimal Physical 
Dipole representing charge distribution of an amino acid that carries a zero net 
charge (C-terminal arginine at physiological pH). The inset image shows the 
plane bisecting the molecule defining the angle plotted in the figure; to illus- 
trate the original charge distribution, the electrostatic potential of the arginine 
is plotted on its molecular surface as a color map with blue representing positive 
and red representing negative potential [8] . The relative error is computed as in 
Fig. 1, (ratio of the absolute error to the RMS average of the reference poten- 
tial), at a distance 2Rq, which as before is twice the distance from the center of 
geometry to the outermost charge. The Optimal Physical Dipole produces less 
than 5% error (solid purple line) for all angles, whereas the optimally placed, 
see [9l , point dipole produces over 30% error (dashed green line) for some 
angles. 

has already been globally minimized. Therefore, to uniquely define the charge, 
q, we follow the 0PM definition and globally minimize the next term in the er- 
ror expansion, namely the octupole term. Specifically, if we consider the £ = 3 
term of Eq. (jl2p . using the connection formula from spherical multipoles to 
Cartesian multipoles we can compute 



where Oijk and Oijk are the octupole moments, in Cartesian coordinates, of 
the original distribution and the Optimal Physical Dipole respectively, for an 
expansion computed about the center of dipole. Now, by noting that Oijk is a 
function of q, we find that Eq. ([2T|) is satisfied when g cxi or if the charge 




(21) 
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value is given by 



3p6 



(22) 



2E. 



Thus, Eqs. dH]), ^ and ^ define the Optimal Physical Dipole, i.e. defines 
the best placement of charges such that RMS error on a sphere centered at 
the center of dipole is minimized. For certain charge distributions relevant to 
molecular biophysics, the Optimal Physical Dipole produces error several times 
lower than that of the corresponding a point dipole, see Fig. [3] 

Although the Optimal Physical Dipole can produce dramatically lower error 
than the point dipole, it is possible that 



In this case, the charge given by Eq. (1^^ is imaginary. This situation occurs 
when the orientation of the dipole with respect to the octupole is such that 
increasing the distance between the charges of the Optimal Physical Dipole 
always increases the error. Thus, Eq. (f21]) is formally satisfied only for q — )■ 
00. In a practical calculation, a physical dipole with the above property is 
constructed by fixing the separation between the charges ||ri — ^2!! to a very 
small value, while increasing the charge accordingly to maintain the original 
dipole moment. In these cases, the Optimal Physical Dipole does not offer 
an advantage over the optimal point dipole, however, the Optimal Physical 
Dipole can always mimic the point dipole to arbitrary precision and thus the 
two distributions will produce equivalent error. Thus, even if inequality (|23p 
holds, the Optimal Physical Dipole represents the optimal placement of two 
point charges and is at least as accurate as the point dipole. 

Curiously, the most biologically important molecule, water, satisfies inequal- 
ity ([23|) . Thus, among point charge representations that use only two charges 
to approximate the true charge density distribution of water molecule, the most 
accurate electrostatic potential outside of the original charge distribution is pro- 
duced by two charges that mimic a point dipole placed at the center of dipole, 
that is itself a special case of the Optimal Physical Dipole with — r2|| — >■ 0. 
Other two-charge charge models that may appear more 'intuitive' by keeping 
the charge-charge distance comparable to the size of the system, i.e. water, are 
bound to produce less accurate representations of the electrostatic potential of 
the original charge distribution. For example, the model used in [12j . places 
a — 0.7|e| charge on the oxygen and +0.7|e| at the geometric center of the two 
hydrogens. However, this 'intuitive' dipole solution produces 20% more error 
at positions near the hydrogen atoms and overall produces ~ 10% higher RMS 
error compared to Optimal Physical Dipole solution (when integrated over a 
sphere of radius 2Ro), see Fig. 01 We stress that the above accuracy analysis 
is strictly limited to electrostatics; other considerations may dictate electro- 
statically suboptimal placements of point charges in coarse-grained models of 
water. 




(23) 
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Figure 4: Relative error in computed electrostatic potential for a water molecule: 
Optimal vs. 'intuitive' Physical Dipole. The reference charge density of the 
electron distribution for water was calculated using the CCSD method with aug- 
cc-pCVTZ basis set [71 HH [23] at experimental equilibrium geometry of water in 
vacuum. The electron distribution is confined to a box with a side length of 4ao; 
the discretization is 0.1 ao along each axis, where ao is the atomic unit of length 
(the Bohr Radius). Here, the Optimal Physical Dipole consists of a negative and 
positive charge with magnitude 10|e| whose center is offset by O.OTOao from the 
oxygen nucleus toward the protons along the water symmetry axis. The charges 
are separated by 0.074ao to maintain the original dipole moment of water. The 
'intuitive' dipole representation, places a positive charge, with magnitude 0.7|e|, 
at the geometric center between the protons and a negative charge (— 0.7|e|) on 
the oxygen. The relative error is computed as in Fig. 1, (ratio of the absolute 
error to the RMS average of the reference potential), at a distance 2Ro which 
in this case has been defined to be 5 times the OH distance from the center 
of the quantum electron charge distribution. Along some directions, e.g. at 
6 = = 2tt, the Optimal Physical Dipole produces error 6 times smaller than 
that of the 'intuitive' dipole. Overall, the average RMS error of the potential due 
to the Optimal Physical Dipole is 10% smaller than that due to the 'intuitive' 
dipole. 



5 Cost Benefit Analysis 

Whether the OPMs introduced in the preceding sections will be a beneficial sim- 
plification in practice will depend upon the specifics of the application. Physical 
multipoles may be easier to implement than point niultipoles in applications that 
already utilize point charges, such as many molecular dynamics codes. Further- 
more, the lowest order Optimal Physical Multipoles have been shown to be 
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always at least as accurate as the corresponding point multipole and in certain 
distributions may be substantially more accurate than an equivalent point mul- 
tipole. We conjecture that the accuracy advantage for some practically useful 
distributions also holds for higher order OPMs. 

One of the major advantages in using the Optimal Physical Multipole method 
versus point multipoles comes in at the level of pairwise interactions. Note that 
although Cartesian formalisms exist for the calculation of pairwise interactions 
[15] , we will restrict the following discussion to the more commonly used spher- 
ical formalism. Consider two systems of charges, one containing of Na charges 
{qi} located at {rj} clustered around an origin located at Ra and another con- 
taining of Nb charges {qj} located at {rj} clustered around an origin located 
at Rb- Then the interaction energy between the two systems is given by 

Na Nb 

UAB=y^y^--^ n (24) 

Thus, the interaction requires NaNb computations. So, for large charge dis- 
tributions the computation scales as 0{N'^). Now, using the point multipole 
expansion, this interaction energy becomes |22| 



oo oo 



2iA + 2lB ^^'^ 



X E E (cX"^(«ab) 



qer^ler^" ^^A'^A'^B,mB\£A + iB,mA + mB) (25) 



for the irregular spherical harmonic, IgJ^'j^ ^^'^ {Rab), given by 

/ . \ 1/2 v-{rnA+mB) ( n \ 

T-mA-mB(T}^ V_ f 47r \ ^iA+lB y^-AB) 

^^^^>-\2tA + 2tB + l) ■ ^^^^ 

where Rab = \\Ra ~ Rb\\, Rab is the unit vector in the direction of Rab 
and the bracketed expression is the Clebsch-Gordan Coefficient. This expres- 
sion has the benefit of scaling linearly with the number of charges. However, 
the number of computations required to compute the pairwise interaction en- 
ergy between two distributions scales as 0{l^) with the number, £ of point 
multipole terms in the interaction. In particular, to determine the interaction 
energy between two charge distributions each with zero net charge but non-zero 
dipole moment requires the calculation of 9 terms, namely the dipole-dipole 
interaction term. On the other hand, an equivalent Optimal Physical Dipole 
representation requires only 4 calculations. To determine the interaction energy 
of a system with non-zero monopole, dipole and quadrupole moments, a scheme 



14 



based on point multipolcs would require, in general, 84 terms, namely the vari- 
ous monopole-monopole, monopole-dipole, monopole-quadrupole, dipole-dipole, 
dipole-quadrupole and quadrupole-quadrupole terms. On the other hand, the 
interaction energy of two Optimal Physical Quadrupoles, containing 5 charges 
each, requires only 25 terms. Additionally, the Coulomb interaction, used to 
compute the interaction between Optimal Physical Multipoles, does not require 
the computation of spherical harmonics or Clebsch-Gordan coefficients, which 
may be cumbersome. Thus, using Optimal Physical Multipoles has a compu- 
tational advantage over spherical point multipoles when calculating interaction 
energies. We have kept the above discussion at the most general level specific 
symmetries present in the original charge distribution may change the numbers. 

6 Conclusion 

Point multipole expansion is a widely used approach to simplify and approx- 
imate potentials and fields around complex charge distributions via power se- 
ries. Hoever, in practical calculations which typically retain only a few leading 
terms, the approximation can lead to inaccurate results in the near-field. In this 
work, we have introduced an alternative to the point multipole expansion for 
approximating the electrostatic potential of a charge distribution, namely. Op- 
timal Physical Multipoles (OPMs). An 0PM consists of a set of point charges 
which are optimally placed to best reproduce the original charge distribution 
at a given miiltiople order. Namely, the 0PM exactly reproduces the point 
multipole expansion of the original distribution up to a given order while opti- 
mally approximating the remaining lowest order non-zero point multipole terms. 
Thus, by construction, OPAIs retain many of the useful properties of point mul- 
tipole expansions, in particular they retain the asymptotic behavior of the point 
multipole expansion. At the same time, an expansion based on OPMs can be 
more accurate than the point multipole expansion of the same order; for some 
charge distributions relevant to molecular biophysics (amin-acid) the accuracy 
gain is more than 5-fold in the near-field. At least for the lowest order OPMs 
(monopole and dipole) , the expansion is guaranteed to be at least as accurate as 
the corresponding point multipole expansion of the same order. We have pro- 
vided a general framework for calculating OPMs to any order and have derived 
closed-form expressions for the Optimal Physical Monopole and Dipole. 

In comparison to point multipoles, expansions based on Optimal Physical 
Multipoles have many desirable properties that may be useful in practical com- 
putations; in particular, their mathematical form - the sum of contribution from 
point sources - is simpler that that of the conventional point multipole expan- 
sion. Thus, OPMs may be easier to implement into existing molecular dynamics 
protocols. At the same time, unlike many approximations based on point charge 
representations, OPMs are uniquely defined and preserve the natural hierarchy 
of the multipole features of the original charge distribution. 

Optimal Physical Multipoles is a new concept; thus its many applications 
and potentially useful properties remain unexplored in this proof-of-concept 
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work. OPMs are expected to have utility in coarse-grained [31 [T3] and multi- 
scale methods [2] , especially in dynamics [1] where analytic expressions and the 
simplicity of the algorithms is key. The Optimal Physical Multipole framework 
provides a systematic way of deriving approximate charge distributions that 
have the potential to be both computationally effective and produce an accurate 
representation of the original electrostatic potential. To further improve the 
representation of the original potential via OPMs, future work may consider 
partitioning the original charge distribution into several domains, and finding 
OPMs for each of them separately, similar to the distributed multipoles approach 
[TUl [201 m] . Further exploration of the mathematical and physical properties of 
OPMs is also desirable; future studies may reveal whether closed form solutions 
exist for higher order OPMs. 
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